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Abstract 

By minimizing free energy density, we show that the stehar core of a hydrogen 
burning star is not in a global thermodynamical equilibrium unless density, temper- 
ature, mass and composition assume given values. The core (as the solar interior) 
may be viewed more appropriately as a metastable state with very long lifetime. 
Slightly non-extensive distribution function could be the natural distribution for a 
weakly non-ideal plasma like a stellar core and represents a more appropriate ap- 
proximation to this system than a Maxwellian distribution, without affecting bulk 
properties of stars. 
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1 Introduction 



Landau and Lifsits [1] distinguish between two different equilibrium states 
called metastable and stable states. A system in a metastable state cannot re- 
turn to its initial state after a sufficiently large deviation from it and could pass 
to another (meta)stable state. Although metastable states are stable within 
certain limits, after a time that could be very long, the system will pass to 
another state. A system removed from the stable state with maximum entropy 
will sooner or later return to it. 

Sewell [2] shows that the equilibrium states are globally thermodynamically 
stable if they minimize the free energy density of the system, satisfy KMS 
(Kubo-Martin-Schwinger) fluctuation-dissipation conditions and are of Gibb- 
sian type. Metastable states (ideal with infinite lifetime or normal with finite 
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long lifetime) have a minimum of free energy density only in some reduced 
space state and these states are always supported only by models with many- 
body, long-range interactions among the particles (as, for instance, gravita- 
tional. Coulomb, spin-glass interactions). Considering infinite systems, Sewell 
minimizes density of free energy to avoid divergences. The minimization of 
this quantity allows us to analyze local, finite subsystems of the entire system 
(like, for instance, as small as the Debye-Hiickel plasma sphere). 

Recently, it has been shown that metastable states can naturally be described 
by distributions accounting for temperature fiuctuations [3,4] or by quasi- 
particle models [5] , and belonging to the generalized statistical mechanics that 
now is usually called non-extensive thermostatistics [6] . 

In this work we are interested in non-relativistic, hot, dense electron-nuclear 
plasmas (like the core of stars or the solar interior) and discuss the conditions 
required for their equilibrium or their metaequilibrium (metastability), follow- 
ing the approach outlined by Landau and Lifsits [7] and recently applied by 
Vasiliev [8], but minimizing free energy density and total energy density. 

It is commonly accepted that stars like the Sun have a core, i.e. an electron- 
nuclear plasma, in equilibrium conditions, that can be described by an electron 
and an ion velocity distribution which are Maxwellian (local thermodynamical 
equilibrium) . 

Main sequence stars produce a luminosity equal to the heat production rate 
and reach a stationary regime. Nuclear fusion reactions produce the heat and 
do not sensibly perturb the state of the Sun. This stationary state that is 
usually called, in astrophysics, equilibrium state, has a long life-time, until the 
nuclear fuel is completely burned out. Stars are not in thermal equilibrium at 
present because they radiate and produce heat in the core with a temperature 
decreasing from core to surface. Thermal exchanges take a very long time, 
whereas gravitational equilibrium is rapidly reached. Nuclear reactions are so 
slow that in the time interval between one fusion reaction and another the 
system returns to the Maxwellian equilibrium due to many binary Coulomb 
collisions among the ions [9] . 

We want to point out that the return to the initial state after a fusion reaction 
is also possible if the distribution function describing the initial state is a 
stationary (stable or metastable) one, analytically expressed by a deformed 
Maxwellian (e.g. a Druj^venstein distribution) or a non-extensive distribution. 
What counts to decide the formal expression of the distribution is the type of 
collisions in which the ions are involved and the dependence on momentum 
of the elastic collisional cross sections (Coulomb, screened Coulomb, enforced 
Coulomb, among others), or the complete absence of ion- ion correlation and 
interaction [10]. Therefore, this argument based on the slowness of fusion 
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reactions is not useful to decide whether the system (the electron-nuclear 
core) is really in an equilibrium state (globally thermodynamically stable) or 
is close to it, or rather it is a metastable state, and to decide what distribution 
is the best approximation. 

In this work, we deduce the values of temperature, density and core radius of 
stellar interiors having a given composition at which they are in equilibrium 
(or in metaequilibrium state), supposing stellar cores made of ideal ionized 
gas plus corrections due to particles identity and internuclear interactions, as 
indicated in Refs. [7,8]. 

We find that globally thermodynamically stable equilibrium is reached in stars 
of the same size and composition of the Sun, at the core thermal energy of 
IcbT" Ri 5 keV (more than three times higher than the actual temperature of 
the solar interior), and at an electron density of n ~ 0.027- 10~^^fm~^ (about 
half the actual value in the solar core, that must have a volume of the central 
core of radius 0.2 Rq). 

Therefore, Maxwellian distributions are of course very good approximations 
(zero approximations), considering the high precision reached by SSM (Stan- 
dard Solar Model) in describing solar neutrinos production. However a bet- 
ter approximation should be the distribution of a metastable state, possibly 
given in an effort to accommodate equilibrium and metaequilibrium states, 
as recently developed in non-extensive thermostatistics, in spite of criticism 
devoted to its use in this field [11,12]. 

As already shown in the recent past, very small deviations from Maxwellian 
momentum distribution do not modify the properties of stellar core, are in 
agreement with the helioseismology constraints, but may affect the evaluation 
of the nuclear fusion rates that may be enhanced or depleted, depending on 
super diffusion or subdiffusion property of the particles [13]. Let us mention 
that deformed distributions have been recently used also in the interpretation 
of many heavy ions processes in thermalization and equilibration of quark- 
gluon plasma [14]; a study of big bang nucleosynthesis of lithium and of a 
possible solution of the problem of lithium abundance in astrophysical systems 
have been undertaken by us and is in progress, among other applications [6]. 



2 Metastable states 

Effort to accommodate a fundamental characterization of either equilibrium 
or metastable states in terms of certain stability conditions was undertaken 
since long time. The equilibrium states are globally thermodynamically stable, 
and are those that minimize the free energy density of the system and satisfy 
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KMS conditions. Grades of metastability correspond to different local stability 
properties. 

Metastable states (ideal) do not minimize free energy density of the entire 
system, but minimize the restriction of free energy density to some reduced 
state space. They are long-lived, supported by long-range, many-body forces. 
They are locally thermodynamically stable if the free energy of the system 
cannot be lowered by modifications of the state that is confined to bounded 
regions of space [2] . 

The metastability of a state depends on the thermal perturbations to which the 
system is subjected. A state may have very long lifetime in presence of certain 
thermal perturbations and may be short-hved when subjected to others. Of 
course, in presence of perturbations or thermal fluctuations, the distribution 
function of the system differs more or less from the Maxwellian distribution. 

By using well known results of Landau and Lifsits [7], and following the ap- 
proach of Vasiliev [8] , we examine the conditions for equilibrium by minimizing 
the free energy density. The free energy of an ideal electron gas is 



where N and m arc, respectively, the number and the mass of electrons, n is 
the electron density, k^T is the thermal energy of the plasma. 

Let us consider now a hot dense electron-nuclear plasma: the free energy should 
be written, starting from Eq. (1), in the following form 



where -Fquantum and -Fcorrciations accouut for the main corrections that arise 
in real gases, namely the quantum exclusion principle between fermions and 
correlation effects between charged particles. 

The first correction, Fquantum, can be written as 




(1) 



ideal ~l~ -^quantum ~l~ -^correlations ) 



(2) 



quantum 



+NEin , 



(3) 



with 
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where ao ~ 0.53 • 10^ fm is the Bohr radius, and ei ~ 2.68 • lO^MeV^/^fm^. 
The second correction is 

i^correlation = -NE2n^^^ , (4) 



with 



J (kBT)V2 - (kBr)V2 ' 



where Z is the average chemical composition and 62 — 2.04 MeV^^'^fm^. 
If we define 



a — 



mc^keT 
27r(;ic)2 



3/2 



ai(kBT)3/2 , 



where ai ~ 0.28 • 10"^ MeV~3/^fm~^ the free energy density, starting from 
Eqs.(2), (3) and (4), is 



F 
V 



-keTn - keTnln (^^^ + Eiu^ - E2n^''^ , 



V being the volume of the plasma. 

Now, we should calculate the value n* of particle density that minimizes the 
F/V free energy density: thus, we can write the extremal condition 



d{F/V) 



dn 



0, 



T,N 



(5) 



from which we obtain the non-linear equation 



In 



n 



ai(kBT)3/2 



+ 2 



ei 



3 (^„v. = o, 



■ n 



(kBT)3/2 2 " (kBr)3/2 



(6) 



that links together T, n and the chemical composition Z. 

It can be shown that condition (5) is not only necessary, but also sufficient to 
assure that the root of Eq.(6) is the actual minimum of free energy density. 
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provided that the inequahty 

^ < 9 17 • 10^ MeV"^/^ 

(kBT)3/2 ^ 



is satisfied. Astrophysical plasmas of stellar cores usually have k^T ~ 1 keV 
and Z of order unity: then the previous inequality holds for these systems. 

Looking for a second relation linking n and T, we can take the total energy 
of non-relativistic electron-nuclear plasma as 

^totai + £;kinetic + ^NE,n - -NE2n}'^ + ^ [ j T/keT , (7) 



where the last term represents the black-body radiation contribution. 

According to virial theorem [15], the potential energy U of particles with 
Coulomb interaction is equal to their double kinetic energy with opposite 
sign, namely 

[/ = -2 ^Nk^T = -SA^keT . (8) 



According to Eqs. (7) and (8), we obtain the total energy density 

^totai _ 3 3 e, _ 3 e.(Z+lp ±{Wrf 

V ~ 2^^2(kBT)V2^ 2 (kBT)V2 ^ ^15(;ic)3 



whose minimum condition. 



0, 



Ny 



can now be written as 

2"" 4(kBr)3/2^ +4 (kBr)3/2 ^ ib\nc) 

Eq. (9) is the second equation that links T, n and Z together. 

Considering that, in a Sun-hke star, n ~ 10-^^ fm"^ kBT - 10-=^MeV and 
Z ~ 1.25, we can neglect the terms with ei and 62 in Eq. (9), thus obtaining 

1/3 

= 6^3^1/3 ^ (10) 



k^T 



3 /fc 
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where b ~ 4.56 ■ 10^ MeV^fm^ 



Therefore, k^T and n must satisfy both Eq. (6) and Eq. (10), given a fixed 
chemical composition, Z. If we choose Z — 1.25, i.e. the Sun composition, we 
find out from the previous discussion that the plasma is globally thermody- 
namically stable if and only if its electron density n* and its thermal energy 
keT* are, respectively, 

~ 2.74 • 10-^^ fm-^ and k^T^ ~ 5 • 10"^ MeV . (11) 



Let us now consider a star characterized by the values of and k^T^ given 
above in Eq. (11); we want to determine if its mass, M*, could satisfy the 
condition ^ Mq, where Mq ~ 1.99-10^^ g is the solar mass. We suppose the 
star in hydrostatic equilibrium, namely the gravitational pressure, which tends 
to compress the star, is exactly counterbalanced by the internal pressure of 
the hot gas and by the pressure of radiation. Therefore, the static-equilibrium 
equation is 

GM' _ , TT ^kBT)^ 



with V (xR^. 

If we impose that n = n^, k^T = keT* and M = M* = M©, we obtain the 
following relation between the radius of a stellar core in equilibrium conditions, 
i?*, and the actual solar radius, Rq, 

R^ ~ O.21i?0 . (12) 



By comparing the quantities n, T and M, we have seen, through Eq. (12), 
that the core of a star like the Sun is not in an equilibrium state globally 
thermodynamically stable, although the equilibrium Maxwellian distribution, 

used for solar core description, is a good approximation. The core can be 
better described as a metaequilibrium or mctastable state. This state may be 
due to temperature fluctuations [4] or quasi-particle behavior of the ions of 
plasma [5]. 



3 The number of peirticles inside the Debye-Hiickel sphere 



Before reporting on the application to main sequence stars (with masses close 
to the solar mass) of free energy density minimization, we wish to recall the 
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fully analytical descriptions of stellar models and stellar stability by Clay- 
ton [16], Nauenberg and Weisskopf [17], and Balian and Blaizot [18]. 



As reported in these papers, the equilibrium properties are obtained by hy- 
drostatic equations with the addition of local thermodynamical equations for 
matter and radiation, based on Boltzmann-Gibbs entropy, with the conse- 
quence that stars evolve along a sequence of states in a complete hydrostatic 
and thermal equilibrium. However, as we shall verify, in these states physical 
quantities like density, temperature and core radius of stars with the com- 
position and size of the Sun are not those required to achieve a global free 
energy density minimum, that, rather, is a local minimum with a very long 
lifetime, because of many-body interactions. Following Sewell, this state is a 
metastable state which (as shown for instance by Sakagami and Taruya [19]) 
can more conveniently be described by a non-extensive distribution function. 

The star interior is a hot dense electron-nuclear plasma at high pressure and 
temperature. Density and temperature are growing toward the center, because 
we require that the system be stable in the gravitational field. At zero approx- 
imation, plasma core can be considered as a Boltzmann ideal gas, with energy 
?>N'k-s,T/2. SSM, based on ideal gas approximation of the core [20,21], cor- 
rectly evaluates the measured solar neutrino fluxes. SSM cannot account for 
the existence of density perturbations since it is based on hydrostatic evolution 
equation. 

In the hydrodynamical approximation, density perturbations can be induced 
by corresponding temperature fluctuations due to convection of matter be- 
tween layers with different local temperatures. 

Greater is the quantity (ksT)^/^?!"^/^ (where n is, as usual, the number den- 
sity), better approximation of the stellar core is the equilibrium Maxwellian 
distribution. In fact, this quantity is related to the square number A^d of par- 
ticle inside the Debye-Hiickel sphere, which must be sufficiently large in order 
to apply Boltzmann-Gibbs statistics. 

If we define the volume of the Debye-Hiickel sphere Vd, 



and assuming that Vb contains only protons, helium nuclei and electrons, for 




and the number of particles A^d inside the Debye-Hiickel sphere. 



A^D = nVo oc 



(kBT)^/^ 

nV2 
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a star with composition similar to that of the solar core, we obtain A^d ~ 4. 
This is a very small number to apply, without criticism, the Boltzmann-Gibbs 
statistics, because the fluctuations of that number are great. Furthermore, the 
value of Ad can increase if k^T increases and/or the density n decreases. 

Going on with time (few billion of years) the core of stars evolves toward states 
that are closer to equilibrium, but they deviate from this behavior when the 
helium burning becomes the most important fusion process. At t ~ 8 • 10^ y 
we have an electron-nuclear (hydrogen) plasma with a density lower than 
the actual density and a temperature much higher [22]; therefore this state 
(pre-helium burning state) is closer to global equilibrium and the Maxwellian 
distribution describing this state is a better approximation than at an earlier 
time because (kBT)^/^n~^/^ is greater. 

A state with a great value of (kBT)^/^n~^/^ is closer to equilibrium than a 
state with a lower value. Let us also note that, although Rj^ oc (keT/n)^/^ is 
constant along the star profile, (kBr)^/^n~^/^ is not. 

Of course, as Njj increases, the weak electron screening correction 

J-WES f ^ ^1^2 \ 



goes to one. 

In non-extensive thermostatistics the free energy of an ideal gas (without 
considering corrections at the moment) is given by [23] 



Arf(g)l^^°(kBT) 



\3/2AfE 



where 



1 f m \3^d/2 
At 



^ Yd! Utt/I^J 



lug X = (a;^~* — !)/(! — g) is the g-logarithm and d{q) is a suitable function of 
the entropic parameter. Let us apply the above expression of Fq to the Debye- 
Hiickcl sphere. When q 1, in order to have a minimum of Fq, one needs 
a large value of Rd that means a large value of (kBT)^/^n~^/^. Of course, as 
5 — > 1 the free energy reaches the minimum for fixed T and n. 

The entropy of the subsystem can be deduced from the relation 
F — TS — E'tot , 
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the expression of F being well-known. 

Neglecting two terms in comparison to others, we can write 



— ~ -iVo 
kB 2 



+1 + n 



ei 



(kBr)3/2 



(kBr)3/2 



(13) 



We know that the non-extensive Tsallis entropy can be written, for a system 
with a small Nd, as [24] 



D 



+ 1 + 



l-g5 
2 2^ 



D 



(14) 



Equalling the two relations in Eqs. (13) and (14) and substituting n* and T* 
with A^D ~ 4, we obtain that in the Debye-Hiickel sphere the q parameter must 
assume the value of ~ 1.001, in order to achieve equilibrium. The correction 
of one thousand, compared to unity, seems very negligible. However it can be 
not at all without importance in the calculation of the fusion rates. 



4 Distribution function of metastable states 



Non-equilibrium systems with metastable states can be composed of regions 
that show fluctuations in space and time of an intensive quantity like the 
inverse temperature /? [4]. 

Non-equilibrium metastable states of a macroscopic system are made of cells 
that are temporarily in local equilibrium. Each small region has a different /3, 
varying in accordance with the probability distribution function f{(3). There- 
fore, stationary probability density of a non-equilibrium system are given by 





The function /(/?) can assume several expressions. However, in case of small 
deviations from Maxwellian distribution, the P{s) function can always be 
written as [4] 

l + ^(g-l)/?oV + ...] , 



P(£) = 



-f3oe ^-&{Poef 



— e 



where 5 — {1 — q) /2, q is the Tsallis non-extensive parameter, and (5q is the 



10 



average value of the quantity p. In this case, if the distribution has a 
standard deviation a, we obtain the useful relation 



P{e) is proportional to non-extensive Tsallis distribution when the entropic 
parameter q is very close to unity. 

For a pure self-gravitating system Sakagami and Taruya [19] have shown that, 
extremizing non-extensive Tsallis entropy, a stellar polytrope distribution can 
be obtained. This polytrope function is regarded as a finite isothermal poly- 
trope, usually as an equihbrium state. However, the one-particle distribution 
function of stellar polytrope shows velocity spatial dispersion. Therefore it is 
no longer the equilibrium but the quasi-equilibrium state (metastable state). 

A stellar core with n, T and M different from the n*, and equilibrium 
values can more properly be described by a non-extensive distribution of a 
metaequilibrium state, rather than by the Maxwellian equilibrium distribution 
(that can be seen as a zero approximation). 

Temperature fluctuations AT/T of order between 1% and 8% have been in- 
voked to analyse neutrino oscillation parameters [25,26]. It is straightforward 
to show that fluctuations of such an order of magnitude are related to values 
of 6 of about 10^"^ 10^'^. With the use of these small deformations, we do 
not change the macroscopic features of the core and its bulk properties; with 
the modified nuclear rates, we evaluate neutrino fiuxes in agreement with the 
measured ones and, within actual uncertainties, these deformations should be 
taken into consideration for a more precise determination of the oscillation 
parameters. 

This kind of fluctuation could be ruled out in an electron-nuclear plasma, after 
analysis of experimental results and global fits. Fluctuations of other intensive 

quantities can be responsible of metastable states. Fluctuating quantities can 
be, for instance, effective chemical potentials or functions of fluctuating energy 
dissipation in the flow (for the turbulence application). 

Other effects that can be responsible of deviations from global thermodynam- 
ical equilibrium should be investigated, as reported below. 

Li and Zhang [27] have discussed how solar core may be far from equilibrium, 
but rather in a quasi-stationary state, investigating the self-generated mag- 
netic flelds in the solar core [28]. In thermal nuclear fusion plasmas, where 
inelastic ions collisions (e.g. pp chain) produce thermal energy, dissipative ef- 
fects of thermal Bremsstrahlung-radiation loss-emission cause a stellar core 
state which differs from a global equihbrium state. Elastic collisions frequency 
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of the order of 10^ Hz would destroy collective motion and coherent structure 
in the plasma if the solar system were isolated adiabatically. In the case of 
an isolated system, the Sun would be in thermodynamical equilibrium. How- 
ever, the solar system is an open system and the nonequilibrium effect cancels 
the collision effect. Consequences are that nonequilibrium can be a source of 
order, the existence of random radiation fields (random electric microfields 
distribution and magnetic field distribution) arc due to self-organization phe- 
nomena. For a discussion of the role of random electric microfield component 
see Refs. [13,29,30]; see also Jiuhn Du [31] for a discussion on self-gravitation 
systems. 

Grandpierre [32] has discussed how the core is neither a ball of ideal gas, nor 
a quiescent steady-state fusion reactor, but a complex self-organizing system. 

Burgess et al. [33], among others, report on the importance of the resonance 
between g-modes and magnetic Alfven waves in the solar radiative zone, how 
strongly density variations affect the solar neutrino survival probability and 
how oscillation parameters depend on the magnitude of solar density fluctua- 
tions. 

The presence of at least one of these effects can influence the statistical distri- 
bution function and be the cause of deviations from Maxwellian distribution. 
We argue that these effects can be taken into account by the appropriate value 
of the entropic parameter q [34] . 



5 Conclusions 

The system we have considered in this work is the electron-nuclear plasma of 
the core of Sun-like stars. We can deduce that the larger is the magnitude of 
the quantity (kBT)^n~^/^, the nearer the system is to the globally thermody- 
namically stable equilibrium state, because the number of particles inside the 
Debye-Hiickel sphere increases and reaches its maximum as 
the maximum value. In this case, the use of the Boltzmann-Gibbs statistics 
is fully justified and the free energy has a minimum; also the free energy of 
non-extensive statistics has a minimum. 

We have looked for the minimum of free energy density as stated by Sewell 
to find the global equilibrium states. Considering the corrections to an ideal 
gas due to identity of particles and to inter-nuclear interaction and the black- 
body radiation emitted, by minimizing the free energy density of the electron- 
nuclear plasma, we have obtained the following values 

~ 2.74 • 10-^^ fm-^ , kei; ~ 5 • 10"=^ MeV and ^ 0.2Rq , 



12 



with a typical chemical composition Z = 1.25. 

States with different values of k^T (lower) and n (higher) are metastable 
states that can be featured by temperature fluctuations or density fluctuations, 
by quasi-particle models or by the presence of self-generated magnetic fields 
or random microfields distributions. These physical effects are responsible of 
metastable distributions that are given by expressions that may be derived 
from the non-extensive thermostatistics. 

Deviations from momentum Maxwellian distribution are quite small for the 
astrophysical systems considered in this paper and the main properties of 
the stars do not change at all for these deviations. Nevertheless, very slight 
deviations can sensibly affect the evaluation of nuclear fusion rates (as we have 
shown elsewhere both for nuclear astrophysics processes and for heavy ion 
physics and quark-gluon plasma properties) and could be useful in solving the 
problem of lithium abundance in the universe, among many other applications. 



Supported in part by MIUR-PRIN2003. 
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